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Abstract 

We reexamine a recently proposed non-inflationary solution to the 
monopole problem, based on the possibility that spontaneously broken 
Grand-Unified symmetries do not get restored at high temperature. We 
go beyond leading order by studying the self-consistent one-loop equations 
of the model. We find large next-to-leading corrections that reverse the 
lowest order results and cause symmetry restoration at high temperature. 
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1 Introduction 



Since the pioneering works by Kirzhnits and Linde fl|], Weinberg || and Dolan and Jackiw 
||, the analysis of the phase diagrams of gauge theories at finite temperature has been 
object of intense research. Many of these studies are connected with spontaneously broken 
gauge theories where the typical picture which emerges is that symmetries get restored 
at high temperatures. 

In standard early universe scenarios based on Grand Unified Theories (GUT's), one 
then starts from the symmetric phase of the model, and as the temperature is lowered, 
several phase transitions occur until the broken phase described by the Standard Model is 
reached. A generic feature of this process is the unavoidable formation of topological de- 
fects during these phase transitions and, as it is well known, the over abundant production 
of monopoles at this stage represents one of the most serious drawbacks of GUT's. 

Different solutions to the monopole problem have been suggested along the years. 
Apart from inflation, which is the most popular one, several authors have proposed alter- 
native scenarios. Among them, let us mention the work of Langacker and Pi 0], who have 
suggested that the existence of an intermediate phase, with no unbroken U(l) symme- 
try, would cause a rapid annihilation of the monopoles produced during the GUT phase 
transition. Another and more radical possibility has been considered by Salomonson, Sk- 
agerstam and Stern || and very recently by Dvali, Melfo and Senjanovic ||. It is based on 
the attractive idea that the monopole producing phase transition could have not occurred 
at all. In other words, it is argued that, no matter how high the temperature was, the 
symmetric phase of 577(5) was never realized. 

Although not very popular and somehow counter-intuitive, the phenomenon of sym- 
metry non-restoration has been present in the literature for quite a while. In fact, it 
dates back to the classic paper of Weinberg |J and some interesting phenomenological 
applications have been discussed in the past 0. 
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The possibility of symmetry non-restoration in GUT's is analyzed in ref. for the 
minimal 577(5) model, which has two Higgs fields in the representations 5 and 24. Ne- 
glecting the gauge couplings, it was shown there that there exists a range of parameters 
for the Higgs potential, leading to the symmetry of the Standard Model at low energies 
but such that the heavy Higgs responsible of the SU (5) breaking keeps a non vanishing 
vev even at very high temperatures. Gauge couplings are considered, instead, in ref.||, 
where it is shown that, due to the large value of the gauge-coupling constant at the scale 
of the 577(5) phase transition, symmetry non-restoration in the minimal model requires 
such large values for the Higgs parameters that perturbation theory breaks down. It is 
subsequently observed that if the Higgs sector is enlarged with a scalar in 45 representa- 
tion of SU(5) symmetry non-restoration can be achieved for values of the Higgs coupling 
constants which seem small enough as to make reliable the lowest order computation on 
which the whole analysis is based. 

In this paper, we will reanalyze the phenomenon of symmetry non-restoration for the 
GUT models considered in ref. |§ by taking into account next-to-leading order corrections. 
The reasons for doing this stem from our study of the simpler O(Ni) x 0(N 2 ) scalar 
model || (which mimics the scalar sector of GUT's) where we found that the inclusion of 
sub-leading corrections reduces in a significant way the region of parameter space where 
symmetry non-restoration occurs. As we will show in this paper, the effect of including 
these corrections is even more dramatic here as they reverse the lowest order result and 
cause symmetry restoration in the entire range of parameters that was considered in ref. || . 

The paper is organized as follows. In Section 2 we present the model and discuss the 
phenomenon of symmetry non restoration to lowest order, while in Section 3 we present 
the results of our next-to-leading calculations. Finally, we leave for Section 4 some final 
remarks. 
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2 A review of symmetry non-restoration 



In this Section, we shall briefly review the basic results about symmetry non-restoration 
in lowest order perturbation theory. For definiteness, we shall consider the GUT model of 
ref. 0, having the usual fermionic content, namely, three generations of fermions (v&^s*, 
(\I / {)io, / = 1, 2, 3 in the 5* and 10 representations respectively and three Higgs multiplets, 
H24 and X45 hi the complex 5, real 24 and complex 45 dimensional representations 
respectively. 

The model is described by the Lagrangian 

L = + ^(D.HUD^H), + (D^y a (D^) a + {D llX )l{D»x)u+ 

-V{H h ® a , X u)+Lp , (1) 

where, in order to simplify the notation, we have called $ , a = 1, ■ ■ • 5 , Hi, i — 1, ..,24 

and Xu,u = 1, ...,45 the independent components of the Higgs fields $5, H 2 i and %45 
respectively and 

F, v = (d.Al - d v Al - gf k A^A k v )T , (2) 

[T\T j ] = x p k T k , Tr(TT J ) = -5 ij , (3) 

(£> M $) a = 8^ a - igA^)A . (4) 

(Dfj,H)i = d,Hi - igAi(T 3 H ) ih H k , (5) 

(D^u = d^Xu - igA^{T 3 x ) uv Xv ■ (6) 

Here, T 3 H , are the SU(5) generators in the 24 and 45 representations, satisfying the 
same commutation relations as in eq.([|) and such that: 

TriT^n) = c H 5« , (7) 

2 Here and in what follows we shall adopt the following indices convention: Greek letters will be used 
for space-time indices, while Latin letters will denote group indices. Group indices will be denoted with 
the initial letters of the alphabet a, b, c... when running from 1 to 5, with the middle letters i, j, k... when 
running from 1 to 24 and with the last letters u,v,w... when running from from 1 to 45 
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Tr(T^) = c x 8v , (8) 

where Ch — 5 and c x = 12. stands for the fermionic Lagrangian, which we shall not 
write down explicitly as it will play a minor role in our discussion. 

As it is well known, the starting SU (5) symmetry of the model can be broken down 
to the low-energy standard model SU(3) C x U{\)em symmetry in two stages. In the first, 
the field H 24 gets a vev (H) that breaks 517(5) to SU(3) C x SU{2) L x U(l) Y . In the 
second, the light fields $5 and/or X24 break this intermediate gauge symmetry down to 
SU(3) C x U(1)em- It is during the first of these phase transitions that monopoles form 
with the Kibble mechanism. As we explained in the introduction, we are studying here 
the possibility that the phase transition leading to the formation of monopoles does not 
occur, namely that the heavy field H 24 keeps a non-vanishing vev at arbitrarily large 
temperatures. 

According to the mechanism discussed by Weinberg 0, this goal can be achieved if 
one lets the field H 2 i interact, with a negative coupling, with the other Higgs fields $ 5 
and/or \2<i- Such couplings are naturally allowed in the model we are considering. If one 
writes the most general SU(5) invariant renormalizable potential [| as: 

V(H h $ a , X u) = --m^HiHi - m|$:$ a - m\ X * u X u + Kix* u Xu) 2 + 

+ ±\ H (H t H t ) 2 + A $ ($:$ a ) 2 - aixlXuXHiHi) + v(H u $ a , X u) , (9) 

it is for example possible to obtain the desired symmetry breaking pattern at zero tem- 
perature with a v {Hi, <3> a , \u) negligible compared to the other terms. Then the condition 
of boundedness from below on V reduces to: 

\ X \ H > a 2 , X x > 0, X H > . (10) 

3 We are using a normalization for the Higgs fields and coupling constants different from that of ref. || . 
The relation between the two is as follows: m} x — 2m'£, m? H — 4m'^, X x — 4A^, Ah = 16A^, a — 4a', 
where the primed constants correspond to those used in [m. 
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The sign of a is thus unconstrained and one can choose it positive such as to give a negative 
interaction between H24 and %45- Obviously one may also consider the analogous case 
where %45 is replaced by $5 in the a dependent interaction term in eq.@[f|. These are 
precisely the two cases considered in ref.@, and on which we will concentrate our attention 
here. 

Let us briefly review how symmetry non-restoration arises in the model. As it is well 
known, the basic tool to study the symmetry behavior of a theory at finite temperature 
is provided by the effective potential ||. As we are only interested here in establishing 
the existence (or not) of a stable symmetric phase at arbitrarily high temperatures (T ^> 
($ a ), (Hi), (Xu)) an d not in the details of the phase transition, we will only need to 
compute the leading terms, of order T 2 , of the effective potential. This turns out to be 
equivalent to evaluating the leading T 2 — contributions to the self-energies of the Higgs 
fields (Debye masses). 

The result for the Debye masses to lowest order in the coupling constants can be 
obtained by adapting Weinberg's formulae to this model and read |§: 

m 2 (T) = £ x (T,p = 0) = T 2 (^ , (11) 

m%(T) = Z H (T, P = 0) = T 2 (1a*(2 + N H ) - + -/^) = ^ , (12) 

where Nh = 24 and N x = 45 are the dimensions of the representations of the Higgs 
fields and D = 24 is the dimension of the group. The different numerical factors between 
eqs.(|TT|) and (|1^) arise from the fact that xa5 is complex while if 24 is real (the analogous 
result for the case of a coupling between if 2 4 and $ 5 can be obtained by simply substituting 
in eqs. flTTD and (0) to N x and c x the corresponding quantities iV$ and c$). 

These formulae exhibit the basic idea behind symmetry non-restoration: for positive 
ct's the interaction if 2 4 — X45 term gives a negative contribution to the thermal masses 
and might make one of them negative (it is easy to see that condition (|10D implies that 
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they cannot be simultaneously negative). Once this happens the corresponding Higgs 
field maintains a vev at high temperatures and symmetry is never restored. 

As explained earlier, we wish to have (H24) 7^ at high T and according to eq.(|l2l) 
this requires: 

^ ^\ 2 + Nh 3 2 Dc H 
m H (T)<0 =► a>XH -^- + -g— . (13) 

The set of coupling constants satisfying this inequality together with the boundedness 
condition, eq.([l(]), thus represents, to lowest order, the region of parameter space in which 
symmetry non-restoration occurs. 

The presence of gauge couplings plays an essential role: as it is evident from eqs.flllf) 



and fll2)) it conspires against symmetry non- restoration. As a consequence, a has to be 
large enough to overcome their contributions to the thermal masses. But, on the other 
side, if a becomes too large, as a consequence of the boundedness condition eq.(|IO|) A x 
or \h are pushed outside the range of applicability of perturbation theory and the whole 
analysis breaks down. In fact, combining eq. (|T3"|) with eq.(|T0|) one gets: 

1 j 2 + N„ 3 2 Peg V 

The r.h.s. has a minimum for 

X " = 392 N h (n" h + 2) • (15) 
which gives a lower bound for A x and a: 

X x>K = ^(^j^-J , (16) 

A . 2 + N H 
a > a = X H — - — • (17) 

Eqs. (p^6|) (pT7|) make it apparent that there are better chances of keeping if 24 broken at 
all temperatures with small values of the coupling constants by coupling it with Higgses 
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belonging to representations of large dimensions. This remark is relevant because of the 
large value of g 2 , tipycally taken to be g 2 ~ 1/4. With this value for g 2 , one has in fact: 

\ H « 0.15 , (18) 

A x ^0.05 , (19) 
a ^ 0.09 . (20) 
Had one coupled instead if 24 to $ 5 , one would have found: 

A$ « 3.9 , (21) 

£ ^ 0.75 , (22) 

While it is more or less clear that for the 5 representation one needs too large values of 
the coupling constants for perturbation theory to be reliable, there seems to be some hope 
with the 45 representation. Nevertheless the reader should be aware that loop corrections 
may contain powers of iV x A x and thus also for the case of a coupling between if 2 4 and X45 
an analysis of higher order corrections results unavoidable. 

3 Next-to-leading order corrections 

In this Section we will compute the next to leading order corrections to the thermal masses 
of the Higgs fields to determine if the results of the lowest order analysis of the previous 
Section are reliable, when the coupling constants are of the order of those of eqs.( |l8l - f20| ). 



At this point, let us mention that some authors |T(J have argued that Weinberg's re- 
sults on symmetry non-restoration are only an artifact of perturbation theory and that 
in reality symmetry is always restored a high temperatures. Their claims are based on 
the observation that non-perturbative calculations (such as 1/N expansion and Gaussian 
Effective Potential) when applied to O(Ni) x 0(N2) global scalar models give always 
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symmetry restoration. Even though non perturbative, the methods used in [jnj are ap- 
proximations, whose range of validity are not clear to us. We will follow a more standard 
route, by assuming that perturbation theory (conveniently improved) is valid and deter- 
mining whether or not the lowest order results remain true after next-to-leading order 
corrections are included. 

Due to the infrared characteristic behavior of field theories at high temperatures, the 
corrections to the self-energies are expected to be of order e 3 / 2 instead of e 2 , where e is 
the largest among A x , A#, ol, g 2 . A simple way of obtaining these corrections is via the 
analysis of the so-called "gap equations" , which correspond to a one-loop truncation of the 
Schwinger-Dyson equations for the self-energy and which are equivalent to a resummation 
of the ring diagrams of the perturbative series. Alternatively, they can also be derived by 
adding and subtracting from the Lagrangian a temperature-dependent mass counter term 
lTJ, which is then determined self-consistently in such a way as to cancel T 2 -divergent 
terms occurring in the self-energies. This kind of improvement of perturbation theory has 
been discussed recently in the context of studies on the Electroweak phase transition |12| . 
We will write the gap equations equations in the symmetric phase and we will identify 
the region of parameter space for which real solutions for the masses can be found as the 
region of symmetry restoration. 

The one-loop gap equations we are considering are schematically represented in fig.(l). 
Notice that we have neglected the Yukawa couplings between the Higgs fields and the 
fermions. We will comment on this point below. The blobs in the internal lines of the 
diagrams shown there represent the complete thermal propagators of the fields. To the 
order we are interested in, they can be approximated with free-like propagators contain- 
ing the thermal masses instead of the zero-temperature ones. For the Higgs fields this 
corresponds to using (in Euclidean space-time): 

= T^m ■ (23) 
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(Xu(-p)Xv(p)) = 2 J UV 2(T , • (24) 
As for the gauge fields A 1 , their propagators are written in the Landau gauge and read: 

p p i + m L {± ) pr 



where 



and 



/ k k s \ 

(Pt)iiu = <V $rs $su ■ (26) 



A' 2 



(PlU = <V - -£r - ( p ^ • ( 2? ) 

(in eq. r and s denote space indices r,s = 1,2,3). Here m^(T) stands for the lon- 
gitudinal thermal mass of the gauge bosons. The standard high-temperature one-loop 
result for mi{T) in an SU(N) gauge theory with Fj chiral fermions, hj real Higgs and E& 
complex Higgs in the representations Ri, Rj, Rk respectively is given by: 



Specializing this formula to our case, we get: 

m£(T) 2 /5 1 1 1 1 1 \ 23 



of — I — c 5 H — Cin H — H — c& H — c Y ) = — q . (29) 
T 2 y V3 2 2 6 3 3 V 3 y V ; 

(Remember that we have three generations of fermions and that C5 = 1/2, cio = 3/2.) 

In our calculation to order e 3//2 this approximation for m^(T) is good enough. 

The gap equations determining m x (T) and m#(T) to order e 3 / 2 can be obtained by 

evaluating the diagrams of fig. (1) at zero external momenta. A simple computation gives 

the high-T result: 

(E x ) uv (p = 0, T) = (E^ + E° 2 + + E*) , (30) 

where 

(m 2 (T)\ 

E- = AX x (l + N x )T 2 h , (31) 



^ = -2aN H T 2 h I ] , (32) 



J^2 



Dc, 



^ = 2g 2 ^T 2 



N v 



(33) 

1 f — f 

^ = 0, (34) 



with a similar result for (E#)j.,(p = 0,T). Here /?.(?/) represents the function: 

I POO 

%2) = 4W dx ( x 2 + y 2y/2( e (* 2 +y 2 ) 1/2 -l) ' (35) 

which, for small values of ?/ has the asymptotic expansion 

h(y 2 ) = —-— x ^- J—y 2 (la^+j --]+... . (36) 

Upon using this expansion in eq.(|30"D and retaining up to the terms linear in 



x x = ^m x (T)/T 2 and xh = y m H{T) /T 2 one gets the following gap equations: 

2 2 g 2 Dc x 1 + N X N H 

* x = v x-^^ 9 -\^r^ + ^H , (37) 

2 2 9 2 Dc H 2 + N H N x 

where u x , vn are the same as in eqs. (|IlT) (|T2"l). (including higher terms in the expansion 
would lead to corrections to the thermal masses of order e 2 log e, but then one would have 
to consider also two-loops contributions to the Schwinger-Dyson equations, which are of 
this order and that we have neglected in our one- loop computation). 

To lowest order eqs.( |37| - j3"8| ) reproduce the results of eqs.(]llT)(|i"2"|) for the thermal 
masses, but the inclusion of the next-to-leading corrections modifies the equations in 
an important way as they introduce a coupling among the thermal masses of the Higgs 
fields. Notice that this coupling, represented by the term proportional to a, comes with 
a positive sign and then works against symmetry non-restoration. 

Now, eqs. ( [37H381) represent a set of two parabolae xh = Ph{x x ), x x = P x (xh) in 
the x x , xh plane. When these parabolae intersect in the upper right part of this plane 
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(remember that xh and x x are defined to be positive), the mass terms of both Higgs fields 
have a positive sign at high temperature and thus symmetry is restored. The condition 
for the existence of such an intersection is: 



P x {x H = 0) < ^ 



(39) 



where x^ is the positive root of the equation Pn(y) — 0. After some trivial algebra, this 
condition is found to be equivalent to 



a < f{g 2 , X x , A H , a) , 



(40) 



where 



f(g 2 , X x , X H , ®) = H X B h 



2iV v 



2iV v 



N H \ 7T 



3 

H — a 

7T 



A 1 + Nx 
x 4vr 



+ 



A- 



47T 



2 Dc x 
AttN x 9 



Hi/2' 



(41) 



The effect of including next-to-leading corrections can be better visualized by fixing 
the coupling constant A x and showing how the region of symmetry non-restoration in 
the (Xh, cx) plane changes in comparison to the lowest order result. To lowest order, the 
values of coupling constants for which symmetry non-restoration occurs are those in the 
region enclosed by the curve C\ which represents the bound (|I~0|): 



C\ : a = \/X x X 



(42) 



and the curve C2 which represents the condition m 2 H (T) = 0, see eq. 



c 2 



. 2 + N H , 3 2 Dc H 



2N^ 



r n h n x 



(43) 



To next-to-leading order, the curve c 2 is shifted to the curve C3 which is obtained by 
solving with respect to a the equation: 



c 3 : « = f(g 2 , X x , X H ,a) 



(44) 
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The choice of A x is not completely free, since it has to be compatible with the bounds 



that come from an analysis of monopole production in charged particles collisions |13| . In 
ref. || it is shown that in order to fulfill this constraint one needs: 

A, > & (45) 

which is a stronger bound than that given in eq. flTBl) . Figure (2), shows the curves ci, C2, C3 
for g 2 = 1/4 and A x = 71g 2 /135. We see that the curve C3 lies entirely above C\, which 
means that for all values of coupling constants compatible with the bound (|T0|) there is 
symmetry restoration. 

For completeness, we have considered also other values of A x and g 2 . Keeping fixed 
g 2 = 1/4, we have found that symmetry restoration occurs for all values of A x consistent 
with the bound ([16]). The results are instead more sensitive to the choice of g 2 . As it is 
clear already from the lowest order result and our previous work on the ungauged model 
||, smaller values of g 2 make symmetry non- restoration easier. As an indication, Figure 
(3) shows the situation when g 2 = 1/16 and A x = 71/135g 2 . We see that even with this 
unrealistically small value of g 2 symmetry non-restoration survives only in a tiny fraction 
of the lowest order region. On the other hand, even this region disappears for larger values 
of A x . 

As we said earlier, in deriving our gap equations we neglected the Yukawa couplings 
between the Higgs fields and the fermions. However, as they always give a positive 
contribution to the Higgs self-energies, they conspire against symmetry non-restoration 
and so can only strengthen our conclusions. 

4 Conclusions 

In this paper we have analyzed the phenomenon of symmetry non-restoration in a Grand 
Unified Model, beyond leading order. We have used linearized gap equations to determine 
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e 3 corrections to the thermal masses of the Higgs fields at high temperature. For realistic 
values of the SU(5) gauge coupling constant we have found that next-to-leading order 
corrections to the thermal masses of the Higgs fields completely overwhelm the lowest 
order result and cause symmetry restoration. 

An interesting question would of course be to determine the critical temperature and 
the nature of the phase transition in this region of the parameter space. To address this 
issue one should compute the finite-temperature effective potential to an order consistent 



with the linearized gap equations. As it is known [|T^, [14]] a two-loop computation in the 
broken phase of the model is needed to perform this task. Notice that also the form of 
the linearized gap equations would change, due to the presence of additional couplings 
and to the fact that the zero-temperature masses of the fields could not be neglected. 

We have limited our analysis to the particular potentials that were considered in 
ref. 0, which are clearly not the most general £77(5) invariant renormalizable potentials. 
To our knowledge, the region of symmetry non-restoration has not been analyzed for the 
general case even to lowest order. In view of our results, which show that the mechanism of 
symmetry non-restoration does not work in the simple models discussed in ||, an analysis 
of the general case appears now to be necessary, before one can claim that this mechanism 
can provide a solution to the monopole problem. This analysis, although conceptually 
simple, becomes tedious due to the large number of free parameters (consider that the 
general potential, even in the case of the model containing only i7 24 and X45 depends on 
16 real independent parameters). Moreover, as our work shows, it will have to take into 
account also the effects of subleading corrections. 

In conclusion, we believe that the issue of whether the monopole problem can be solved 
in realistic GUT models via the mechanism of symmetry non-restoration remains an open 
and certainly interesting problem. 
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FIGURE CAPTIONS 



Figure 1. One-loop gap equations for the Higgs' self-energies. Solid lines correspond 
to X45 5 dashed lines to H 2 4 and wavy lines to gauge fields. The blobs represent insertions 
of full thermal propagators. 



Figure 2. Plots of C\ (the solid parabola), c 2 (the solid straight line) and c 3 (the dashed 
line), for g 2 = 1/4 and A x = 71# 2 /135. 



Figure 3. Plots of c\ (the solid parabola), c 2 (the solid straight line) and c 3 (the dashed 
line), for g 2 = 1/16 and A x = 71# 2 /135. 
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